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■ Abstract. The properties of Lagrangians affine in velocities are analyzed in a geometric 
way. These systems are necessarily singular and exhibit, in general, gauge invariance. The 

■ analysis of constraint functions and gauge symmetry leads us to a complete classification of 
such Lagrangians. 
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! The Dirac's method of dealing with constrained systems [1-3], developed when looking 

for a way of quantizing such systems by means of a "canonical quantization" -like procedure, 
has been shown to be a very efficient method and the geometric approach to the theory 
of systems defined by singular Lagrangians has deserved a big amount of attention during 
the past years. 

In particular, Lagrangians that are affine in time-derivatives have been analysed [4] in 
the framework of presymplectic geometry and the geometric theory of time-independent 
singular Lagrangians [5]. They are important because their Euler-Lagrange equations 
^ become systems of first-order differential equations, appearing as constraints, instead of 

systems of second-order as it would happen with regular Lagrangians. So they will play a 
relevant role in many cases, not only in Physics, where many equations are first-order, as 
in Dirac equation, but also in other fields as in biology dynamics and Chemistry. On the 
other side, these systems are singular and then they are giving rise to gauge ambiguities 
and gauge symmetries. 

The geometric study of a particular type of singular Lagrangian, those which are affine 
in velocities, was carried out in [4] in the framework of autonomous systems, with the 
aim of studying the inverse problem of Lagrangian mechanics and the theory of non-point 
symmetries from a new geometric perspective. Almost simultaneously Faddeev and Jackiw 
developed a method for the quantization of such singular systems which became soon 
very popular and received much attention from most part of theoretical physicists. This 
procedure of dealing with such systems is usually referred to as Faddeev-Jackiw method 
of quantization [6-12]. 

One of the most important properties of singular Lagrangians is the existence of in- 
finitesimal gauge symmetries which are related to the Second Noether theorem. This, 
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which is particularly important for field theories, has only sense in the framework of 
time-dependent systems, the time playing the role of base coordinates in field theory. 
A geometric approach to Second Noether Theorem was given in [13], but the geometric 
theory of the time-dependent description of such affine in velocities Lagrangian systems 
has never been developed, as far as we know, even if it is very important as the only way 
of fully-understanding the meaning of Noether 's Second Theorem establishing the rela- 
tionship between singular Lagrangians and gauge-transformations. Then we feel that a 
re-examination of the problem of singular Lagrangian that are first-order in velocities will 
be very useful and will allow us a better understanding of the theory. 

The two fundamental aspects of these systems described by these first-order in velocities 
Lagrangians are the presence of first order equations and constraints, which make possible 
this alternative Faddeev-Jackiw method of quantization [6,14]. This method is based on 
the reduction theory for the presymplectic form defined by the singular Lagrangian, but 
it admits an alternative by means of the addition of the constraints with some Lagrange 
multipliers in order to obtain the symplectic extension [11,12,15,16]. Once the symplectic 
structure has been got and, therefore, Poisson brackets have been defined, we can make use 
of the canonical quantization procedure. In particular, the FJ-method uses the reduction 
technique and has been applied in many different fields, running form condensed matter 
[17-19] and astrophysics [20], until fluid dynamics [21] and, of course, high-energy physics 
[22]. Even the own Schrodinger equation can be derived using this method [23]. 

It is also to be noted that the FJ-method can be generalized to include also fermionic 
degrees of freedom, i.e. non-commutative variables [24] and the corresponding canonical 
quantization can also be used in super-symmetric theories [8,25,26] with applications in 
super-gravity (see e.g. [27]). 

The paper is then organized as follows. In Section 2 we summarize the results of [4]. A 
Lagrangian approach to Hamilton equations from a geometric perspective will be given in 
Section 3. The framework for the geometric description of time-dependent singular systems 
as given in [28] will be indicated in Section 4, where we will also include a recipe obtained 
from [13] (see also [28]) for the search of gauge symmetries for singular Lagrangian systems. 
The geometric theory of time-dependent Lagrangians which are affine in the velocities will 
be developed in Section 5 and the gauge invariance of such systems will be studied in 
Section 6. Finally, the theory is illustrated with several examples. 

As a matter of notation, tangent and cotangent bundles will be denoted by tm '■ 
T M — > M and ttm '■ T* M — > M, respectively. The set of vector fields along a map 
/ : M — > N (see e.g. [29]), also called /""-derivations in [30], i.e. maps X : M — > 
T N such that tn ° X = /, will be denoted £(/). Examples of such kind of vector 
fields are T/oF and Z o /, where Y G £(M), T/:TM^TiVis the tangent map 
corresponding to /, and Z G X(N). It has been shown in Pidello and Tulczyjew [30] 
that a vector field X along / determines two /*-derivations of scalar forms on M: one 
of type and degree —1, denoted ix, and other of type d+, denoted dx, defined in the 
following way: given a G /\ P (N) and v\, . . . ,v p -i G T^ M, we define ixot G /\ P_1 (M) by 
(i x a(x), (ui, . . . , Vp-i)) = (a(f(x)), (X(m), f* x M, ... , /**(u p -i))), and d x a G /\ P {M) 
is given by dxOL — dixct + ixda [31]. 
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2. Geometric approach to time-independent 
Lagrangians which are affine in the velocities 

In the geometric description of a time-independent Lagrangian system the states are 
described by points of the tangent bundle tq : TQ — > Q of the configuration space Q, 
which is assumed to be a n-dimensional differentiable manifold. We are interested in the 
geometric study of systems described by a Lagrangian including only terms up to first 
order in velocities, namely, with a coordinate expression 

L(q,v) = m j (q)v j -V(q) , (2.1) 

where summation on repeated indices is understood. 

The property of a function being linear in the fibre coordinates (velocities) of the 
tangent bundle is intrinsic because it is preserved under point transformations. If G 
C°°(T Q) takes the form <f> = 4>i{q) v l in a particular set of local coordinates, then, under 
the change of coordinates q' 1 = q' l (q), we have that in order to 

<M<zK = <^') J 

the components <j>i(q) should transform as 

MQ) = tiW)(^f) > 3 = h..-,n, (2.2) 

i.e., the functions 0j((?)'s should transform as the components of an associated basic 1- 
form a G /\ 1 (T<5), ol = (f>i(q) dq l . This suggests us that there is thus a one-to-one linear 
correspondence between basic 1-forms and linear homogeneous functions that we establish 
next: if fj, G /^(Q) is a 1-form, then ju denotes the function fi, G C°°(T(5) defined by 
ju = itA*, where t G 3£(tq) is the identity map id : T Q — >■ T Q viewed as a vector field along 
the tangent bundle projection tq; in fibred coordinates, t = v l d/dq l . That is, jj,(q,v) = 
(fj,(q),t(q,v)), the contraction making sense because t(q, v) G T q Q. In coordinates, when 
[i = rrij(q) dqi , p,(q, v) = rrii(q) v l . 

The geometric theory for systems described by a Lagrangian L = ft, — TqV , with 
V G C°°(Q), that in coordinates of the bundle induced from those of a chart in Q becomes 
L = rrij(q) — V(q), has been studied in [4]: the energy El and the presymplectic form ujl 
are given by E L = Tq(V), Ol = Tq^-, and, therefore, ujl = —r^dfi), which in coordinates 
reads as follows: 

/ drm dm k \ i k 

u '^{w-w) dq Adq ■ 

The Hessian matrix W with elements Wij = d 2 L/dv l dv^ vanishes identically and, there- 
fore, all the TQ-vertical vectors, i.e. ^ l {q,v) d/dv l , are in the kernel of ujl- 

The search for the other elements in the kernel of ujl starts by looking for a local basis 
{Z a = (zl)}, of the module of eigenvectors (if any) corresponding to the null eigenvalue of 
the matrix A with elements given by elements 



(2.3) 
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Then, a basis for the kernel of ujl is generated by the vector fields X a = z l a (q) d/dq l and 
d/dv l , with a = 1, . . . ,n — no, where no denotes the rank of the matrix A. The primary 
constraint sub manifold M is then determined by the constraint functions (see e.g. [5]) 
4>x = XEl, with X E KeruL, that in the case we are considering are 4> a = X a E^, i.e. 

<Pa = zl(q)^, (2.4) 

because the energy is TQ-projectable, El = TqV ', and then dE^/dv 1 = 0. 

The general solution for the dynamical equation i(X) ujl = dE^ will be given by 

X =W + * a 4]^ + f^- i , (2-5) 
with rj l (q) being a solution of 

dV 

An rf = g- , (2.6) 

and where A a and f % are arbitrary functions on T Q. 

There is a special class of vector fields D in T Q which are called second order vector 
fields, hereafter shortened as SODE fields, which are characterized by S(D) = A, where S 
denotes the vertical endomorphism [32-34] and A is the Liouville vector field generating 
dilations along the fibres. They can also be characterized by Ttq o D = t. 

The constraint functions for a Lagrangian given by (2.1) are basic functions, i.e. holo- 
nomic constraints, defining a submanifold Q' of Q. Consequently, the secondary constraints 
functions for the existence of a solution restriction of second order vector field will be d4> a 
and are given by linear functions in the velocities. A solution X given by (2.5) can be the 
restriction of a SODE only in those points of T Q for which 

dV 

In these points, the general solution of the dynamical equation is given by 

while the SODE condition corresponds to the choice A a = 0. 

The particularly simple case in which [Q,dn) is a symplectic manifold, i.e. det A ^ 0, 
and therefore Q is even dimensional, was also studied in [4], where it was shown that 
then Kera;^ is made of all TQ-vertical vectors and therefore there will be no dynamical 
constraint function. The globally defined solution of the dynamical equation is 

^ dq 1 ^ dv 1 
with the functions rf uniquely determined by 



(2.9) 
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The Marsden-Weinstein theory of reduction for the presymplectic system defined by 
the Lagrangian (2.1), (TQ,ul = — Tg(d/x), El = Tq(V)), establishes that the reduced 
symplectic manifold is (Q, —dfi) with Hamiltonian function V. 

The Hamiltonian formalism for the Lagrangian (2.1) was also studied in [4]: the primary 
constraint submanifold Pi = JFL(TQ) is determined by $>j(q,p) = pj — rrij(q) = 0, i.e. is 
the graph of the form fx. Using the identification of it with the base Q the pull-back of the 
canonical 1-form ujq in T* Q is — d\i. The Poisson brackets of the constraint functions are 
{$-,-, = {pj — m,j(q),pk — rn>k(q)} = Ajk, and therefore, when d\x is symplectic, all the 
constraints are of the second class. The Hamiltonian function H is defined on Pi by the 
restriction of the function V = KqV. The general theory leads again to the study of the 
Hamiltonian dynamical system (Q, — dfx, V), as in the Lagrangian case. 

Let now L be a singular Lagrangian for which all constraint functions O = XqEl, 
a = 1, . . . ,n — no = N, are holonomic; we consider an extended configuration space 
WL N x Q and denote by pr 2 the natural projection pr 2 : M. N x Q — > Q. We can then 
introduce a new Lagrangian Li in T(M Ar x Q) by 

Li = L + A a </> a , (2.10) 



where the tilde stands for the T pr 2 -pull-back and A a are the new additional coordinates 
(whose corresponding velocities will be represented by (, a ). Taking into account that 

w Ll = 2i, E u =E L -\ a 4> a , (2.11) 



we see that Kerc^Li is generated by the set of vector fields projecting onto the vector fields 
X a of KeruL, plus d/d\ a and d/d( a . The constraint functions determined by d/d\ a are 
the (pull-back of the) original primary constraint functions 4> a and SODE condition leads 
us to consider the tangent bundle of the new configuration space M fc x Q' . The solutions 
of the dynamics will project under pr 2 onto the solutions of the original problem. 

A similar approach can be done when we use d<p a instead of <p a as constraint functions 

and we replace the original Lagrangian L for L2 = L + X a d(j) a . In this case E-^ 2 = El and 
c^l 2 = ujl + dtfia A d\ a , from which we can see that 

^ 9 .„ d ,9 „ d 
x = Cin + f +'/V? + i - 



dq { dX a dv i d( a 
is in Ker cul 2 if an d only if 



M - £ a V0 a - Wti = 0, e-V0 a = o, wz = o, 

from which we see that we will obtain as constraint functions (the pull-back of) those 
obtained directly from L, and the dynamics will correspond, up to the gauge ambiguity in 
the coordinates A a , to the dynamics obtained in Q'. 

The relation X a d(f> a — d(X a 4> a ) — ( a 4> a shows that the Lagrangian L 2 may be replaced 
by L 3 = L — C a O 5 which is quite similar to Li with the change of A a for its velocity ( a . 

In the more general situation for which d/i is singular, the primary constraint functions 
(2.4) will be holonomic and therefore the previous remarks show us that they will determine 
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a submanifold Q' C Q characterized by some constraint functions <p a and the corresponding 
secondary constraint functions for a second order evolution will be d<f) a . Every such a 
constraint can be incorporated in a new Lagrangian L defined in the tangent bundle of a 
new configuration space Q of the form Q = Qx IR( n ~ n °), by 

L(q\ A, v\ = L(q\ v*) + X a d^ a (q\ v% (2.12) 

where (q l , A a , v l , ( a ) denote the coordinates on the tangent bundle T Q. 
The expressions for 6>l and cul are 

#l = 9l + \ a d0 a = ujl + d<fia A d\ a , (2.13) 

and therefore the rank of wl may be higher than that of ujl and this is the starting point 
in the Faddeev-Jackiw approach. The energy is the pull-back of E L: = E L . 

3. A Lagrangian approach to Hamilton equations 

Let M be the configuration space of a mechanical system and consider Q = T* M 
endowed with its exact canonical symplectic structure u> = —dOo, where 6q is the Liouville 
1-form in T* M (see e.g. [35]). Then, given a function H e C°°(T* M), let us define the 
linear Lagrangian L e C°°(T(T* M)) by 

L = 9q—t^* m H, (3-1) 
which in local coordinates is written 

L(q,P;v,u) =pjV J - H(q,p). (3.2) 

In this case 9 L e C°°(T(T*M)) is given by 9 L = r-f. M 9 and E L = t^, m H. The 
matrix A given by (2.3) is now the symplectic matrix 



J = 



/ 
-I 



The Kernel of ujl = —d9L is then made up by the tt* M-vertical vector fields 

Y f, g = /*(?. p; «) ^? + 9 % (q, p; w ' M ) ' ( 3 - 3 ) 

and therefore the presymplectic system defined by L, (T(T* M), ujl, El), admits a global 
dynamics which is not uniquely defined but given by 

OH d dH d , 9 , d 

where / and (7 are arbitrary functions. The integral curves of each one of such vector fields 
are determined by the system of differential equations 

( .i dH 

dpi 
dH 

^ % dq i 
u l = f*(q,p;v,u) 
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However, only in those points for which 

.- dH 



* = % (3 ' 5) 

dq 1 

the solution can be chosen to be the restriction of a SODE vector field. The preceding 
equations determine a submanifold C on T(T* M) and the condition on the restriction of 
the vector field r f i9 to be tangent to C determines the functions f l and g l by means of 

obtaining in this way the vector field on C 

dH_d__dH_d_ r fdH\_d_ r (^£\JL 

The dimension of C is only twice that of M and then it can be parametrized by (q l ,Pi)- 
The expression of Tc in these new coordinates is 

dH d dH d 



dpi dq 1 dq 1 dpi 

and therefore, here the Hamilton equations arise as determining the integral curves of the 
uniquely defined vector field in the submanifold C in which such a SODE solution of the 
dynamical equation can be found. Then, a curve 7 : I — > M whose lift to T M lies in C is 
the solution we were looking for. 



4. Geometric description of time-dependent singular system 

For the reader convenience we introduce in this Section the notation to be used and a 
summary of several properties and results of interest for the following sections. 

The evolution space of a time-dependent mechanical system whose configuration space 
is the n-dimensional manifold Q is R x T Q [36] , which is the space of 1-jets of the trivial 
bundle 7r: E x Q — > R, j l n = R x TQ; K. is endowed with a volume form dt e f\ (K) 
and represents the Newtonian time. There is one vector field in R, d/dt G X(R), such 
that i(d/dt)dt = 1. The main geometric tools to be used in the geometric description 
of time- dependent mechanics are those of jet bundle geometry [37,38]. The fc-jet bundle 
of 7T is J k n = R x T k Q, with T k Q representing the space of k- velocities. In particular, 
J 2 n = R x T 2 Q is the space of accelerations and J 1 n = MxTQ, the space of velocities. 
By convention we will write J° = R x Q. For each pair of indices k, I such that k > /, 
there is a natural projection Tik,i : R x T fc Q — >■ M, x T l Q and we will denote iTk = 
7r o ii k : R x T k Q — > R, the projection of J k n onto R. 

If a e Sec(7r), then j k a G Sec(7Tfc) will denote the fc-jet prolongation of a. So, if 
a(t) = (t, o-*(t)), we have j k a(t) = (t, a l (t), da l /dt, . . . , d k a l /dt k ). We also recall that the 
differential 1-forms 6 G f\ 1 (J k ii) such that (j k a)*6 = 0, no matter the section a, are 
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called contact 1-forms of J tt. They are the constraint 1-forms for the so called Cartan 
distribution. 

The theory of time-dependent Lagrangian systems makes an extensive use of the notion 
of vector field along a map. In particular, there exist vector fields along irk+i,k, T fc , 
representing the total derivative operators. They are defined by means of T k o j k+1 a = 
T(j k a) o d/dt, V<7 E Sec(Tr). For every F E C°°(R x T h Q), the Lie derivative of F with 
respect to T fc , i.e. the function d T kF = i Tk dF E C°°(R x T fc+1 Q), represents the total 
time derivative of F, usually written as dF/dt or simply F. In fibred coordinates (t, q l , v l ) 
for R x TQ and the corresponding ones, (t, q l , v l , a 1 ), for R x T 2 Q, we have T° = d/dt + 
v i d/dq i E X(7Ti 5 o) and T 1 = d/dt + v i d/dq l + a} d/dv l E X(tt 2i1 ). Obviously T° and the 
operator t E 3c(tq) introduced in Section 2 are related, T p o T° = t o p 2 , with p and pi 
being the projections onto the second factor of R x Q and Ix TQ, respectively. 

A vector field X E X(R x Q) can be lifted to R x TQ giving rise to a unique vector 
field X 1 E X(R x TQ) which is 7ri -projectable onto X and preserves the set of contact 
1-forms of E x TQ. If the coordinate expression of X is X = r d/dt + X 1 d/dq l , then 

Vector fields of type X 1 are called infinitesimal contact transformations (hereafter ict); it 
is worthy to note that 

(fX) 1 = (7Tl f)X 1 + fX V , 

where X v = S(X 1 ). Here S is the vertical endomorphism on Ix TQ [31], which is a (1,1)- 
tensor field whose expression in the fibred coordinates (t, q l , v l ) is S = d/dv l <8> (dq l — v l dt). 
Note that the local 1-forms given by 9 l = dq l — v l dt generate the set of contact 1-forms of 
R x T Q and the Cartan distribution is but Ker S. 

A similar definition works for the prolongation of vector fields X = r(t, q,v) d/dt + 
X l (t, q, v) d/dq l E X(7r lj0 ): the vector field along n 2 ,i given by 

is the first prolongation of X. Thus T 1 = (T ) 1 . 

More details about these notions and constructions can be found in [28]. 

The key object on which the geometric formulation of the dynamics corresponding to 
a time-dependent Lagrangian L E C°° (RxTQ) is based is the Poincare-Cartan 1-form, 
defined by 

6 L = dLoS + Ldt E /\\r x TQ). (4.1) 

In fibred coordinates Ol = (dL/dv l ) 9 l + Ldt. Another important object related with L is 
the Euler-Lagrange 1-form, which is defined by SL = i T idQL E /\ (R x T 2 Q), with local 
coordinate expression SL = Li9 l , where 

_ d (dL\_ (dL\ 
Ll ~ dt{dv*J 2,1 W/ 

are the variational derivatives of the Lagrangian L. 



On singular Lagrangians affine in velocities 



9 



The dynamical equation to be considered here is 

i(T)dQ L = 0, (4.2) 

and, according to Hamilton's principle, the condition for the section a G Sec(7r) to be an 
extremal of the action is 

(j 1 a)*[i(Z)dQ L } = , VZGl(lxTQ). 

The problem is to find a vector field V G £(R x T Q) which is solution of the dynamical 
equation (4.2) and whose integral curves are the first prolongation j 1 a of sections a G 
Sec(7r). In other words, V must be a time-dependent SODE, i.e. a vector field T such that 
T 7Ti 5 o ° T = T°. Its integral curves are parametrized by t and the local expression of V is 
T = d/dt + v l d/dq i + P d/dv\ Both sODE-type and iCT-type fields belong to Ker S, but 
in general a SODE is not an ict. 

The submanifold of R x T Q where the dynamical equation (4.2) possesses such kind 
of solutions is given by the following theorem [28]: 

Theorem. Let Ml denote the coisotropic subbundle 

M L = {Ze T(R x TQ) | S{Z) G Vi l0 (Kerde L )} , 

and Vi 5 o( - ) means the ixi^-vertical part, which is but the kernel of the differential of the 
Legendre transformation TL: IxTQ — > R x T* Q, that is, Vi ; o(Ker c/©l) = KerTjFL. 

Then, there exist solutions of the dynamical equation which are restrictions of a SODE 
field in the points, and only in that points, of the set defined by 

M = {zeRxTQ \ dQ L (D,Z)(z) =0, VZg M l , D any sode }. 



When L is regular, Ml = {0} and M. = R x TQ and, consequently, there is no 
restriction on the motion. But there exist (primary) constraints for a singular Lagrangian 
which are given by the following conditions 

<$> z = dO L (D,Z) = 0, Z G Ml- (4.3) 

The functions $^ are the primary constraint functions. Obvious conditions for the consis- 
tency of the dynamics V compatible with the constraints (4.3) are given by 

Xz = T($ z ) \ M =0 (4.4) 

which either give rise to the secondary constraints or (partially) fix the dynamics V. When 
the process is iterated, we will hopefully arrive at the final constraint submanifold -M/, on 
which there exist solutions of the dynamical equation which are the restriction onto M. / 
of SODE fields tangent to Mf. 

Note that here the only ingredient is the singular Lagragian L which provides both 
the (nonholonomic) constraint functions and the dynamics. In this sense this is a problem 
that does not coincide with the more frequently studied constrained situation in which the 
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starting point is a given nonholonomic constraint distribution [39]. This latter situation, 
that is receiving much attention during the last years (see e.g. [40,41]), and is important 
by its relation with the theory of connections and by its control theoretical applications, 
is different from the present case of first order singular Lagrangians. 

The dynamics obtained by applying the constraint algorithm sketched above may be 
not unique, a fact which is related with the gauge invariance of the Lagrangian. The 
appropriate geometric tool to deal with a gauge infinitesimal transformation 

where g = g(t) is an arbitrary function of the time, is that of a vector field along the 
bundle map nip. In fact, let {X a = X l a (t, q, v) djdq 1 \ a = 0, 1, . . . , R} be a family of 
R + 1 7r- vertical vector fields along 7Ti o and g(t) an arbitrary function in R. Then, if X g 
is the 7r- vertical vector field along 7Ti i0 

X o = e E > ( 4 - 6 ) 

its first prolongation Xg is the infinitesimal generator of the gauge transformation (4.5). 
Such X g is said to be a gauge symmetry of the Lagrangian L if there exists a function 
F g EC°°(RxTQ) such that 

d xl L = XlL=^. (4.7) 

In such case (5L 7 X g ) + dG g /dt = 0, with G g = F g — (Ql,X 9 ), and conversely, if there 
exists a function G g such that (5L 7 X g ) + dG g /dt = 0, then d x iL = dF g /dt with F g = 
G g + (0_l,X 3 ). (The contractions (5L,X g ) and {®L,X g ) make sense because of the 7r 2) o- 
and 7Ti o-semibasic character of 5L and G^, respectively.) 

All this mathematical apparatus can be used to give a geometric version of Noether's 
Second Theorem as it can be seen in [13,28]. The theorem essentially establishes that a 
gauge-invariant Lagrangian is necessarily singular and it satisfies the so called Noether 
identities 

El -,\a da ( SL i X a) _ n 
1 J dt a 

a 

i.e. 

E, ■,\a d a (LiX l a ) _ 
[ ' dt* 

with Li being the variational derivatives of the Lagrangian L. However, a singular La- 
grangian needs not to be gauge-invariant and a method for the determination of the gauge 
symmetry underlying a given singular Lagrangian has been developed in [13] (see also 
[28]). The method is based on the necessary conditions which are derived from the sec- 
ond Noether's theorem and it also tells us whether, or not, a given Lagrangian is gauge- 
invariant. Let X g given by (4.6) be the wanted gauge symmetry of L and choose a vector 
field X a e X(R xTQ) such that T7r lj0 ° X a = X a . Assume that 

R-l 



G 9 = J2 G 



d a g 
a dF 

a=0 
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It follows from the symmetry condition (4.7) that the vector field Xr belongs to the 
distribution Ml — Ker S and the functions G a satisfy the relations 

Gi? = 0, G a + G Q _i + (5L,X Q ) = (a = l,...,R), G + (SL, X ) = , 

from which it follows that all the functions G a must be jFL-projectable and satisfy the 
recursive relations 

G Q _! = -de L {D,X a ) - D(G a ) , 

with D being any SODE; in particular, Gr-i = —cIQl(D, Xr) is a primary constraint 
function (see (4.3)). The algorithm for the determination of the gauge symmetries proceeds 
by determining in an iterative and orderly way the functions G a and the vector fields X a 
along the following steps: 

1. Choose Gr = and select X R E (M L — Ker S) fl Ker T tt\ in such a way that Gr-i = 
—d<S>L(D,Xfi), with D being any SODE, is a ^"L-projectable primary constraint func- 
tion. 

2. Then, let us determine a ^-vertical vector field Xr-i in such a way that the 1-form 
Xr-i = i(X R-i)d<d l — dG r-i be 7Ti,o-semibasic; Gr-2 is defined by Gr-2 = Xr-i(D). 

3. When the successive functions G a are ^-"L-projectable the process may be iterated and 
when we find Gr-n = Gr = the algorithm enters into a cycle and the solutions 
appear cyclically repeated. We can take R = N — 1 to be the higher order for the 
derivatives of g(t) and the algorithm ends up. 

In the case when in some step there is not any solution, we stop the process and return 
to make (if possible) a new choice for the solution in a previous step. If there is no solution 
in any case, we have to conclude that the Lagrangian is not gauge invariant. 

The next sections are devoted to show the application of these constructions to the 
case of a Lagrangian which is linear, or more accurately affine, in the velocities v l . 

5. Geometric theory of time-dependent 
Lagrangians which are affine in velocities 

From the geometric viewpoint, a time-dependent Lagrangian L which is affine in the 
velocities arises from a 1-form A G /\ 1 (1R x Q) in the following way: L = i T o\ e C 00 ^ x 
T Q). In fact, if A = TOj(£, q) dq l — V(t, q) dt, then the Lagrangian, to be denoted by A, is 

X = m i (t,q)v i -V(t,q) . (5.1) 

Obviously the time-independent case we have dealt with in Sect. 2 is simply a special case 
of this: given \i E /^(Q) and V E C°°(Q) then we consider the 1-form A = p* \x — p*V dt E 
/\ 1 (M x Q) which yields the (p 2 -pull-back of the) time-independent Lagrangian (2.1). 
Coming back to the general case, the basic geometric features for this Lagrangian are: 

1. The Poincare-Cartan 1-form is given by 6^ = tt^qX. 

2. jFA-projectability means yr^o-projectability, because 

FX = fx o A o 71-1,0 , 
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where ji = ttr x id-p* q is the natural projection of T*(R x Q) onto R x T* Q. 

3. Vi, (Kerrfe^) = KerTyr^o, so that = T(R x TQ). 

Taking into account all these facts, we find that the primary constraint functions are 
given by 

<&z = dQ^(D, Z) , ZgT(IxTQ), DanySODE, (5.2) 

i.e. $z = i{Z){i^odX). The primary constraint manifold is described in the following 
terms: 

M = {zeRxTQ \ <$> z (z) =0} = {zeRxTQ \ i T od\(z) = 0} . (5.3) 

However, if we recall that there exists a local basis of £(R x TQ) is made up from a sode 
and vector fields Y 1 and Y v = ^(y 1 ), with Y G Ker T n [28], we see that the only effective 
constraint functions are those given by 

<5>y = dO^Y 1 , D) , Y G Ker T n . (5.4) 

Then the functions $y are in a one to one correspondence with the elements of Ker T it. 
They also verify the property 

$ Yl+fY2 = $ Yl + KV)^, / g C°°(R x Q) . (5.5) 

We will say that $y and $y are 'linearly dependent' if <E>y = $jy for some / G C°°(Rx Q) 
everywhere non-null. This property trivially takes place when Y = fY but in the case 
when Y and Y are not dependent it is a property related with the existence of a gauge 
symmetry of the Lagrangian as we will see later. 

In local coordinates (t, q l ) for R x Q and (t, q\ v l ) for R x T Q, respectively, the vector 
field Y is written Y = Y l (t, q) d/dq l and dX = \ Aij dq l A dqi — Ui dq l A dt, where 

_d mj dim _dV drrij 

Aj ~ d? dqi Ut ~ d? + dt ■ (5 - 6) 

Consequently, the primary constraint functions are 

$ y = {A ij v i -Ui)Y\ (5.7) 

i.e. they are affine in the velocities. A basis for such constraint functions is made up from 
the following functions 

$i = ^d/dqi = Aij >v J - Ui = 0, i = 1, . (5.8) 

Note that these equations are the Euler-Lagrange equations obtained from the Lagrangian 
(5.1) and all of them appear in this formalism as constraint equations. 

We can deduce from (5.4) that a primary constraint function <Ev is 7ri ^-project able 
(i.e., holonomic) if and only if the 1-form i(Y) dX is 7r-semibasic. In an equivalent way, $y 
is 7Ti )0 -projectable if and only if 

Y G fl\ = [Ker T n] _LdA n Ker T n . 
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(The superscript _L dX means <iA-orthogonal complement). In local coordinates as above, 
fi> is spanned by vector fields Y = Y l (t, q) d/dq x such that A^Y 1 = and then the 
corresponding constraint function $y = u>i Y l = is projectable. 

Obviously, the maximum number of linearly independent constraint functions equals 
the dimension n of Q and the maximum number of the holonomic ones is p = n — rank A, 
where A is the skew-symmetric matrix of elements A^ (5.6). 

As far as the dynamics is concerned, any SODE T is a solution on M. of the dynamical 
equation because of (5.2). The consistency conditions (4.4), r$y = 0, will generate addi- 
tional (secondary) constraints and/or fix (maybe partially) the dynamics. More accurately, 
the consistency conditions for non-projectable primary constraint functions will determine 
some components of the SODE, while the projectable ones will give rise to secondary con- 
straints that are affine in the velocities, which fix the dynamics. The uniqueness is got 
when all the constraints so obtained are independent and in the maximum number. The 
two fundamental cases are. 

I) The first fundamental case arises when fi A = 0. Then the matrix A is regular and 
therefore there are no primary holonomic constraints. The n independent constraints (5.8) 
are non-holonomic and lead to n equations determining in a unique way the dynamics on 
M, i.e. 

v j = (A-y^k . (5.9) 

Of course, this is only possible when the dimension of Q is even, n = 2m. 

Let us analyse this "regular" case in geometrical terms. The condition Q\ = means 
that the distribution Ker dX is 1-dimensional, so dX is a contact form on RxQ generating an 
exact contact structure onlxQ. In fact, every nonzero vector field Z G Ker dX satisfies the 
condition i(Z)dt ^ 0, as it follows from the trivial fact that Z also belongs to [Ker T7r]- LdA , 
and if the condition i(Z)dt = is fulfilled Z should also be 7r-vertical and, consequently, 
Z = 0. Moreover for every pair Z, Z' G Ker dX the vector field Z" = [i(Z)dt] Z'-[i(Z')dt] Z 
lies in fi A , Z" G fi A , i.e. Z' = [i(Z')dt/i(Z)dt]Z . 

On the other hand, no new constraints arise from the consistency conditions (4.4) and 
the n independent non-holonomic constraints will fix the dynamics on M. in a unique way. 
Let r\ G KerdA be a vector field such that i{rf)dt = 1. It generates locally Ker dX and in 
coordinates turns out to be 

„.! + (,!-)««,£. (5.10) 

The description (5.3) of M. means that T°(z) G Ker dX (TTifi(z)), z G M., and, consequently, 
we can assert that the ICT t] 1 is, at least on M, a SODE field. It is such that i(r] 1 )d<3)^\M = 0, 
thus the dynamical equation (4.2) is satisfied by i] 1 on M and the SODE T compatible with 
the constraints is the one generated by i] 1 . 

In summary, when Q\ = 0, the Lagrangian system on R x T Q reduces to the (in 
general, time-dependent) Hamiltonian system (R x Q, dX, H). The manifold R x Q is the 
extended phase-space and the Hamiltonian function H is essentially the energy E^, which 
is a holonomic function: £~ = nl H; in local coordinates H = V. 

Both dynamical systems are equivalent in the sense that the integral curves of the 
dynamics in R x T Q are the first prolongation to R x T Q of the integral curves of the 
Hamiltonian dynamics in R x Q. 
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Using the 2-form dX we can define a Poisson bracket in 1R x Q according to the following 
construction. Given a 1-form « G ^(1 x Q) there exists a vector field X a e X(K x Q) 
such that z(X a )<iA = a iff i(Z)a = 0, VZ e Ker <iA; obviously such a X a is not unique, the 
indeterminacy being Ker dX itself. Although a does not satisfies the condition above one 
can take the "<iA-semibasic" part, given by 

a x = a - ^ dt, Z e Ker dX , (5.11) 
%[Z) dt 

which, in fact, depends only on Ker<iA and not on a particular Z e KerdA, and it trivially 
annihilates KerdA, i.e. there exist vector fields X a such that i(X a )dX = a x . 
In local coordinates, if a = ciq dt + ai dq\ then 

a x = cii (dq l - (A- 1 )^ ojj dt) 

and 

X ce = X° V +(A- 1 y k a J ° 



dq k 

where X° is an arbitrary function and r\ is given by (5.10); in particular, for a = df, 

d x f = (df) x = ^w-(A- 1 y^ J dt) 

and 

The Poisson bracket {/, g} of the functions / and g is then defined by the rule 

{f,g} = dX{X g ,X f ) = d x g(X f ) . (5.12) 

In local coordinates, 

the fundamental Poisson brackets being {g J , q k } = (A -1 )^ (see [42]). So the equation 
of motion derived from the constraint equations (5.8) can be written in the form qi = 
{q^q k }u k . 

In the autonomous case (Q, —dy) is a symplectic manifold and we obtain the Hamil- 
tonian system (Q, — dfx, V) analysed in detail in Section 2. In particular, the equations of 
motion read 

dV 

(see (2.9)). 

II) Second case: Q\ ^ 0. Then A is singular and there will be r = rank A indepen- 
dent non-projectable constraint functions, and then the number of effective projectable 
constraint functions is not greater than p = n — r. These functions generate secondary 
constraint functions that together with the non-projectable functions will determine, at 
least partially, the dynamics. The uniqueness of solution for the dynamics depends on 
whether or not the primary and secondary constraints are independent. The analysis of 
gauge invariance in the following section will clarify these points. 
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6. Gauge in variance of first order Lagrangian systems 

In this section we will show how to make use of the algorithm of gauge symmetry 
explained in Section 4. It starts by taking Gr = and choosing a vector field Xr G 
(M^— Ker 5 , )flKer T tti in such a way that the function Gr-i given by Gr-i = cI<S)^(Xr, D) 
be a TTi^-projectable primary constraint function. Then, we can choose Xr = Y 1 (that is, 
Xr = T TTi o ^" 1 = ^°7Ti,0) where Y G O^), so that Gr-i turns out to be Gr-i = $y that 
is a 7Ti o-projectable function. Moreover, we can see that the 1-form A^_i = i(XR-i)d<d^ — 
cIGr-i is 7ri o-semibasic, no matter of the choice of Xr_i. We can write Xr-i as a sum 

X R _, = Y 1 1 +Y 2 v , y 1} y 2 eKerTTr; (6.1) 

and, consequently, 

Gr-2 = d&z(Y?, D) - D{<$> Y ) = ® Yl - X Y (6.2) 

is the difference between a primary constraint function and a secondary one. The algorithm 
only works if it is possible to choose Gr-2 as being a 7ri ^-project able function. The 
iteration of this procedure will give rise to a sequence of functions GR-k analogous to 
Gr-2, and therefore a gauge symmetry will be obtained in this way if we arrive at a 
secondary constraint function xy that it was a primary one, making possible to have 

G R -2 = 0. 

Thus the starting point to have a gauge symmetry is to know the set of holonomic 
(primary) constraint functions. If two (or more) of them are linearly dependent it will be 
possible to choose Xr in such a way that Gr-i = 0. In the opposite case we need to 
know whether or not the primary constraint function gives rise to a secondary constraint 
function which is a primary one. 

More specifically, there are the following four possibilities and correspondingly the four 
types of affine in velocities Lagrangians: 

I) Q\ = 0, i.e. A is regular. In this case there are no holonomic constraints and there is 
no gauge symmetry at all. The dynamics is uniquely determined by the non-holonomic 
constraints (see Sect. 5). 

11. 1) Q\ 7^ 0, i.e. A is singular, and there are two (non trivially) linearly dependent 
holonomic primary constraint functions $y and <fry, i.e. $y = <fr/y = (71^ f)$>y I n this 

case, y — fY is non- null and choosing Xr as Xr = (Y — fY) 1 , we will get Gr-i = (i. 
e., the higher order for the derivatives of the arbitrary function g(t) is R = 0) and come 
to the gauge symmetry 

X g = g(Y-fY)o7r 1 , . (6.3) 

11. 2) 0\ 7^ and there is a holonomic primary constraint function $y, 7^ Y G Q\, giving 
rise to a secondary constraint function xy = D(<&y) which is a primary one, namely, there 
exists a vector field Y G Ker Tn such that xy = &y = ^\0^ ■> D)- In this case, choosing 

Y 1 = Y we find that G R - 2 = dQ^Y 1 , D) - D{$> Y ) = 0. Therefore the algorithm tells us 
that R = 1 and the gauge symmetry is 



X g = 9 (Y TTl.o) + 9 ( Y 7Tl,o) 



(6.4) 
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II. 3) Q\ 7^ and none of the secondary constraints is primary. Gr-2 is not holonomic and 
the algorithm cannot go on, i.e. there is no gauge symmetry. The dynamics is uniquely 
determined by the full set of constraints (both primary and secondary). 



In summary, we will have gauge symmetry when there exist holonomic constraint func- 
tions generating secondary constraint functions that generate a free set with the non- 
projectable constraint functions. 

7. Examples 

Finally, several examples will be used to illustrate the Lagrangian analysis made in 
Sections 5 and 6. As a matter of notation we will use subindices instead of upperindices 
in the coordinates q and the velocities v. 

Example 1. The well-known two-dimensional Lotka-Volterra system can be derived 
from the following Lagrangian which is affine in the velocities [43] : 

t ln y lnx < i . m ^ 

L = —v x - —v y - (a lny + blnx -x - y) , 

where a and b are positive constants. Considered as a Lagrangian function onKxT R 2 ^ 
it derives from the 1-form 

A = — dx - — dy - (a lny + blnx -x - y) dt E A* (R x R 2 ,) . 
2x 2y 1 N 



Therefore, 





/n _±\ 



xy 

1 

\xy / 



and u 




The matrix A is regular and then fl\ = 0. All of the primary constraint functions (5.8) 
are non-holonomic, 

^ v v b v x a 

®x = ®d/dx = - - + 1, 3\/ = ®d/dy = + 1, 

xy x 1 xy y 

and, consequently, the reduced system (5.10) onKx R\ is 

The system of differential equations for its integral curves constitutes the 2-dimensional 
Lotka-Volterra system 

x = x(a- y), y = -y(b - x) . 
This is a Hamiltonian system with a symplectic structure 

1 

a = — dy A dx , 
xy 
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i.e. with defining Poisson bracket {x, y} = xy, and Hamiltonian function 

H = a In y + b In x — x — y . 

This Lagrangian is of Type I. 

Example 2. Another regular case is that provided by the Lagrangian studied in [44] 

1 1 

L = g [(^ 2 + 93 W - QlV2 + (?4 - Qi)v3 - 93^4] - ^ ( 2 9293 + 93 + 94) 

coming from the 1-form A EE x IR 4 ) given by 

A = \ [(92 + Q3)dqi - qi dq 2 + (q 4 - qi)dq 3 - q 3 dq 4 ] - ^(2q 2 q 3 + ql + q\) dt , 
from which we obtain 





f° 


-1 


-1 






( ° \ 


A = 


1 










and lo = 


93 




1 










92 + 93 




Vo 





1 






V 94 / 



The matrix A has rank 4 and then, as Q\ = 0, there is not any holonomic constraint. The 
generating set of primary constraints (5.8) is 

$1 = -V2 ~ v 3 = , $ 2 = Vl -q 3 =0, $3 = ^1-^4-92-93 = 0, $4 = ^3-^4 = 0, 

and the secondary ones, T($j) = 0, with i = 1, . . . , 4, determine a unique dynamics which 
is the restriction of a SODE on the constraint manifold M, namely 

_,. _, d d d d d d d d d 

T = T\ M = 777 + 937; 94 7; ^ 94 7^ 92 7^ V V 3 ~ V 4 - h V4- V 2 



dt dqi dq 2 dq 3 dq 4 
The reduced system on R x I 4 is 



' dvi 



dv-2 



' dv? 



dv 4 



d d d d d Tr 

V= w: +93^ 94^ h94^ 92^— e KerdA 

dt dqi dq 2 oq 3 dq 4 

Note that the restriction of V onto M. coincides with that of r\ x . 

Example 3. Let us now consider the Lagrangian of the type II. 1 

L = 91 v 2 + q 2 v 3 + q 2 v 4 - 92 (94 - 93) 

generated by the 1-form 

A = 91 dq 2 + 92 dq 3 + 92 dq 4 - q 2 (q 4 - q 3 ) dt E f\ (Rx IR 4 ), 
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from which we obtain 



A = 



( 10 0^ 
-10 11 
0-100 
V -10 0/ 



and u) 



( \ 

94 ~ 93 
92 



The primary constraint functions (5.8) are 

$i = v 2 , $2 = -vi + v 3 + v 4 - q 4 + q 3 , $ 3 = -v 2 + 92 , $4 = -f 2 - 92 • 

The distribution Q\ is spanned by the two vector fields d/dq\ + d/dq 3 and d/dqi + d/dq 4 , 
and yields to two holonomic primary constraint functions 

$i + $ 3 = q 2 , $i + $ 4 = -92 

that are linearly dependent. The Lagrangian is gauge invariant, and the symmetry vector 
(6.3) and the associated function F g are (R = 0) 



The dynamics V is determined by the set of (primary and secondary) constraints $j, Xii 
i = l,...,4: 

9 d d d d d d 

r = ^T + (f3+t'4-94+93)^— +^3^— +"47 h (C 3 + C 4 - f 4 + V 3 ) 7 hC 3 - hC 4 -— , 

at agi ag 3 oq 4 ov\ ov 3 ov 4 

where C 3 and C4 are arbitrary functions. 

Example 4. The Lagrangian defined by the 1-form 

A 1 

A = (92 - 93) dq-L - 92 dq 3 - (q 2 - 91) 93 dt G /\ (1 
is of Type II. 2. The matrices A and u; are given, respectively, by 



5 



d 



x 






-1 






1 -53 \ 


1 







|. -1 


93 


-1 


1 






V 92 - 91 / 



A = 

The set of primary constraints (5.8) is generated by 

$1 = 2 + f 3 + 93 , $2 = v 1 - v 3 - q 3 , $3 = -v 1 + v 2 - 92 + 9i • 

The rank of A is 2 and therefore the distribution Q\ is 1-dimensional. It is generated by the 
vector field d/dqi + d/dq 2 + d/dq 3 and there is one primary holonomic constraint function, 
namely, $1 + $2 + $3 = 91 — 92, whose corresponding secondary constraint function, 
X = Vi — v 2 , is likewise a primary constraint: x = ®d/dq 1 +d/dq 2 = $1 + $2- Consequently, 
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this Lagrangian is gauge invariant: starting from Xr = d/dqi +d/dq2 +d/dq 3 and choosing 
= d/dqi + d/dq 2 we obtain the gauge symmetry (6.4) 



d d d 

+ 77" + 



dqi dq 2 dq 3 J 



and the corresponding function F g = g(t)(qi — q 3 ) — g(t)q 3 . The dynamics on M. is given 
by 

r = It + ( u 3 + Qs) t^- + (v 3 + q 3 ) -9- + v 3 7^- + (C + v 3 ) t^- + (C + v 3 ) + C- ° 



dt dqi dq2 dq 3 dv\ dv 2 dv 3 

with C an arbitrary function. 

Example 5. The 1-form 

A = t q 2 dq 1 - qi dq 2 - [ 9l - (t + l)q 2 ] q 3 dte f\ 1 (R x R 3 ) 

gives the time-dependent Lagrangian of the type II. 2 

\ = tq 2 v 1 -q 1 v 2 - [qi - (t + l)q 2 ] q 3 . 

The matrices A and u are given, respectively, by 

/ 92 + 93 

oo = -{t+l)q 3 
\ qi -(t+l)q 2 

Two constraint functions of those determined by (5.8) are non-holonomic 
$1 = -{t + l)v 2 ~ qi ~ 93 , $2 = (t + l)vt + (t + l)q 3 , 
while $3 is holonomic, 

$3 = -qi + (t + 1)92 • 

It gives rise to a secondary constraint, X3 = ~~ v i + (t + l)t>2 + 92, which is a primary 
one, %3 = — $i — $> 2 /(t + 1). The constraint function %3 corresponds to the vector field 
— 3/ dqi — (t + l)d/dq 2 and, consequently, A is gauge invariant. The symmetry vector X g 
(6.4) and the function F g are (R = 1) 





The local expression for the dynamical vector field is 

d d q 2 + q 3 d d d 2v 2 + v 3 d d 

r = 777 - 93 q 7 , : q h V 3 q V 3 q 7 , 1 q h C - 



<9t t + 1 ^<?3 <9^i t + 1 dv 2 dv 3 

with C being an arbitrary function. 
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Example 6. A Lagrangian of the type II. 3 is the one provided by a slight modification 
of the 1-form A in Example 4: 

A = (g 2 - 93) <ki + 92 dq s + (qi - q 2 ) 9s dt e^Ix K 3 ) . 

In this case, 

A = 






-1 






( -* \ 


1 







|. -| 




-1 


-1 






\ 92 - ?1 / 



The primary constraint functions given by (5.8) are 

$1 = -^2 + f3 + 93 , $2 = Vi + v 3 - q 3 , $3 = - v 2 + 9i - 92 • 

The distribution Q\ is spanned by the vector field d/dqi — d/dq 2 — djdqz- There is a 
holonomic constraint function, $1 — $2 — ^3 = 2<j»3 — q\ +92, giving rise to a secondary 
one, x = 2^3 — v\ + v 2 , which is not primary, that is, there is no 7r-vertical vector field Y 
such that x = Consequently, there is no gauge symmetry at all, and the dynamics V 
on the constraint manifold M. is unique: 

„ d ( d d \ ( d d 



dt V dqi dq 2 ) V dv\ dv^ 
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